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I. INTRODUCTION 



As is well-known, the radiative B — > X s ^ decay and the semileptonic B — > X s l + l~ 
decay play an important role in the precision test of the standard model (SM) and the 
search for the new physics beyond the SM. Although no evidence of the new physics have 
been found now in experiments, one can put strong constraints on the parameter space 
of various new physics model from currently available data. 

The minimal supersymmetric standard model (MSSM) is the general and most eco- 
nomical low energy supersymmetric extension of the SM. In order to find the possible 
signals or hints of new physics from the date, various scenarios of the MSSM are proposed 
by imposing different constrains on it [l]. Here, we use the pMSSM (phenomenological 
MSSM) model [l[ in our studies. 

When calculating the new physics contributions to the flavor changing neutral current 
(FCNC) processes induced by the loop diagrams involving the new particles, one needs 
some kinds of model-independent parametrization of the FCNC SUSY contributions. The 
mass insertion approximation (MIA) [2| is the best one of such kind of parametrization 
methods. In the MIA, one chooses a basis for the fermion and sfermion states where all 
the couplings of these particles to neutral gaugino fields are flavor diagonal, and leaves all 
the sources of flavor violation inside the off-diagonal terms {A q AB ) i - of the sfermion mass 
matrix, ie., 

(Ml) i . = m\ j + (Al B ) ij (1) 

where m is the averaged squark mass, A,B = (L, R), q = u,d, and the generation index 
i,j = (1,2,3). 

As long as the off-diagonal terms is much smaller than the averaged squark mass, the 
sfermion propagators can be expanded as a series in terms of {5 q AB )^ = (A AB ) i:j / m 2 . 
Under the condition of (A AB ) <C m 2 , one can just take the first term of this expansion 
and translate the relevant experimental measurements into upper bounds on these <5's. 

According the helicity of the fermion partner, the squark mixings can be classified into 
left- or right-handed (L or R) pieces: 

, . > • ( 2 ) 

The LL and RR mixings represent the chirality conserving transitions in the left- and 
right-handed squars, while the LR and RL mixings refer to the chirality flipping tran- 
sitions. Up to now, many interesting works have been done to draw constraints on the 
parameter <5's from the known data. The strong constraints on the mass insertion LL 
and RR, for example, are obtained from the measured AM^, AM^ and ex H, while the 
mass insertions in the LR and RL are constrained by the data of Br(b — > sj) 0,0], the 
ratio e'/e and the electric dipole moment (EDM) of the neutron, electron and mercury 
atomQ. 

Very recently, the possible SUSY corrections to the branching ratios and the CP- 
violatinff asymmetries of B — > Kit, Krj' and <pK decays have been studied, for example, 
in Refs. |4|, la, |6[ . Here, the size and phase of the MIA parameters (S AB ) 3 j i3 with i,j = (1,2) 
play an impotent role in explaining those observed "puzzles" in B experiments or not. 

The measured values of Br(B — ► Xg'j) and Br(B — > X s l + l~) agree perfectly with 
the SM predictions ja,l2|, which leads to strong constraints on various new physics models 
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[lfll . From the well measured B — > X s 7 decays, the magnitude-but not the sign- of the 
Wilson coefficient CV 7 (m&) is strongly constrained. The sign of C 7l {rn b ) and its absolute 
value, however, can be determined from the precision measurements of the B — > X s l + l~ 
decays. The latest Belle and BaBar measurements of the inclusive i? — > X s l + l~ branching 
ratios indicated that the sign of CV 7 (mj) should be the same as the Cf 7 M (m 6 ) H3- There- 
fore, when calculating the upper bounds of the mass insertions (Sj^ B ) 3 j i3 in the MSSM, 
one should not only consider the constraint coming from the branching ratio of B — > X s 7, 
but also the new information for the sign of CV 7 (mfc) from the new measurements of 
Br(B -> X s l+l~). 

This paper is organized as follows. In SeclHJ we show the formulaes needed to calculate 
the branching ratio of B — > Xgj decay at next-to-leading (NLO) order. The new physics 
contributions to the relevant Wilson coefficients are also given in this section. Then 
in Sec JIIIl by comparing the theoretical predictions with the data of B — > Xgj and 
B — > X s l + l~ decay we update the upper bounds on those mass insertion parameters. 
We considered the cases of one mass insertion and two mass insertion approximation. 
Summery is given in the last section. 



II. Br(B X sl ) IN THE SM AND MSSM 

A. Br(B -> X s j) in the SM 

In the SM, the effective Hamiltonian for b — > 57 at scale \i ~ 0(m b ) reads H3 



G F 

^eff = --J=V tb V* s 
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Y^Ci(ti)Oi(ti) + C 77 (m)0 77 (/x) + C 8g (n)0 8g (n) 



i=l 



+ h.c, (3) 



where VtbV t * is the products of elements of the Cabbibo-Kabayashi-Maskawa (CKM) quark 
mixing matrix |l3j. The definitions and the explicit expressions of the operators Oi 
(i — 1 ~ 6, 77, 8g) and the corresponding Wilson coefficients C{ can be found in Ref . . 
In the SM, the Wilson coefficients appeared in Eq. © are currently known at next-to- 
leading order (NLO) and can be found easily in Ref. |l2| . 

At the lower energy scale fi b ~ 0(rrib), the Wilson coefficients at NLO level can be 
formally decomposed as follows 

C^ h )=cf\^) + ^C\ 1 \^ (4) 

where Cf^ and Cf' stand for the LO and NLO order part, respectively. Finally, the 
branching ratio Br(B — ► Xg'y), conventionally normalized to the semileptonic branching 
ratio Br cxp (5 -> X c eu) = (10.64 ± 0.23)%, is given by Q, [3 

Br(B ^ X s7 ) = Br^(B - X c e^)^|^ • [|D| 2 + A + A] , (5) 

|V cb | 2 vrg(z)k(z) 
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with 



DM = c£V 



7 7 



8 



=1 



A(p b 



An 



4tt 

c (o) 



r i( z ) lo g — 



m b 



A(/i 6 



7T 



5^ 

2 

mi 



i,j=l,i<3 

(M&) +ArRe 



cv 0>e// (^)' 



4} 



C 1 °' e// (M6) 



(6) 
(7) 

(8) 
(9) 

where z = (mf le /m p b ° le f, X 1 = 0.5 GeV 2 , A 2 = (m|. - m|)/4 = 0.12 GeV 2 and m 6 /2 < 
A*b < 2m b . The explicit expressions for the Wilson coefficients Cf \ Cj , the anomalous 
dimension matrix 7, together with the functions g(z), k(z), Ti{z) and fij(S), can be found 
in Refs. Q 0, 2|- The term A(fi b ) describes the correction from the bremsstrahlung 



with 



6 1 



Ai 
2 



-A, 
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process b — > S7g, while the term A includes the non-perturbative l/m^ |l7| and 1/ 
corrections. 
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TABLE I: The input parameters entering the calculation of Br(B — » Xg^y), their central values 
and errors. All masses are in unit of GeV. We use Gp = 1.1664 x 10 _5 GeV. 



M w M z 


m t 


m b 


a- 1 




a s (M z ) 


80.42 91.188 


175 ±5 


4.8 


137.036 


0.118 


BR S l m c /m b 


Mb 


sin 2 6\y 


Ai 




A 2 


0.1064 0.29 ±0.02 4.8l£| 


0.23 




-0.5 




0.12 


A A 


P 


fj 




v th v ts 
v cb 




0.854 0.2196 


0.20 ±0.09 0.33 ±0.06 


0.97 ±0.01 





Using above formulae and the input parameters as given in Table HJ we find the SM 
prediction for the branching ratio BriB — ► Xgj), 

Br NLO (B -> X s7 ) = (3.53 ± 0.30) x 10" 4 (10) 

where the main theoretical errors come from the uncertainties of the input parameters 
and have been added in quadrature. By using the same input parameters, it is easy to 
find the numerical values of C 77 (m fe ), V(m b ), A{m b ) and A(m&) as defined in Eqs.fJBJIH]) 

C 7l (m b ) = -0.3052, V(m b ) = -0.0257 - 0.0156/, 
A(m b ) = 0.0033, A(m b ) = -0.0010 (11) 

One can see that it is the Wilson coefficient C-j^(m b ) who determines the branching ratio 
of B — > Xsj decay, and other three terms are indeed very small. 



4 



B. Br(B X a >y) in the MSSM 

At the leading order, the one-loop diagrams involving internal line SUSY particles 
provide new physics contributions to the Wilson coefficients CV 7 , C$ g , as well as the 
Wilson coefficients CV 7 and Cs g of the new operators 0^ 1 and Os g , which have the opposite 
chirality with and 0$ g appeared in Eq.fJHJ). In the SM, the contributions from chiral- 
ffipped operators Oj.y^ are very small since they are strongly suppressed by a ratio 
0(m s /mb)- In the MSSM, however, the contributions from the operators and Os g 
may be not small in the SUSY models with non- universal A-terms 0, 0, 0] • 

At the one-loop level, there are four kinds of SUSY contributions to the b — > s transi- 
tion, depending on the virtual particles running in the penguin diagrams: 

(i) the charged Higgs boson H ± and up quarks u, c, t; 

(ii) the charginos Xi 2 and the up squarks u, c, t; 

(iii) the neutralinos X? 2 3 4 and the down squarks d, s, b; 

(iv) the gluinos g and the down squarks d, s, b. 

In general, the Wilson coefficients at fi\y ~ M\y after the inclusion of various contri- 
butions can be expressed as 

CM = + Cf + cf + Cf° + Cf, (12) 

where Cf M , , Cf + , Cf and Cf denote the Wilson coefficients induced by the penguin 
diagrams with the exchanges of the gauge boson W ± , the charged Higgs H ± 1 the chargino 
X121 t ne neutralino Xi 2 34 and the gluino g, respectively 1 

In principle, the neutrolino exchange diagrams involve the same mass insertions as 
the gluino ones, but they are strongly suppressed compared with the latter by roughly 
a ratio of a/a s ~ 0.06. The charged Higgs contribution are proportional to the Yukawa 
couplings of light quarks and relevant only for a very small charged Higgs mass and very 
large tan/3. Therefore, We will concentrate on the chargino and gluino contributions only. 
The SUSY contributions to the Wilson coefficients CV 7i 8 9 and CV 7i 8 9 induced by various 
Feynman diagrams at the my/ scale have been calculated and collected, for example, in 

Refs. ESElIIII. 

First, we consider the chargino contributions. As emphasized in Ref. 0, the chargino 
contributions induced by magnetic-penguin, and chromomagnetic-penguin diagrams de- 
pend on the up sector mass insertions (^l)32,3i and (^x,) 32,31, while the LR and RR 
contributions are suppressed by A 2 or A 3 , where A = 0.22 is the Cabibbo mixing. 

As for the mass spectrum of the squarks, the authors of Ref. p considered two cases: 
(a) all squarks have the same mass m; and (b) the stop-right tit is lighter than other 
squarks. But the numerical results show that the difference between these two cases are 
rather small. We therefor consider the first case only. 



1 The parameter At = VtbV t * in Eq.(10) of Ref.[6( has been absorbed into the definition of Cf and Cf 
here. 
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At the high energy scale fiw ~ Mw, the chargino part of the Wilson coefficients, 
C^mw) and Cg g (mw), are of the form p 

C^(m w ) = X; 1 {[(8l L ) 32 + \{SIM Rf + m L ) 32 + X(5 U RL U Y t R* L ] , (13) 
Cl{m w ) = AT 1 {[(S u LL h2 + X(5 U LL U Rf + [(^) 32 + A(^ L ) 3 i] YtRf } , (14) 

with 



Xi), 



R i L 9 = El^il 2 ''xwiPffai) -Y b J2VnU i2 x Wi ^P* 

i=l i=l ^ 

2 

R ™ = -J2 v * v £ x ™ p , L t( x i), ( 15 ) 



8=1 



and 



x (-2 - 9x + 18a; 2 - 7a; 3 + 3x{x 2 - 3) logfz]) 
9(1 -x) 5 ' 



LR x (13 - 20x + 7x 2 + (6 + Ax - Ax 2 ) log[x]) 

7 I 3 ') ~~ 



6 1- x 



.7' 



-1 + 9x + 9x 2 - 17x 3 + 6x 2 (3 + x) log[x]) 



3 v 7 12(1 -x) 5 



.r i 



■1 - Ax + 5x 2 - 2x(2 + x) log[x]) 



9 K ' 2(1 -x 4 



where Y t = \f2m t /{vsm.(3) and Yj, = \f2mb\/l + tan 2 (3/v are the Yukawa coupling of 
the top and bottom quark, xy/ l = m w/ m Xi i Xi = m x t /™' 2 ' ^ = ™ 2 / m Xi with i = (1,2). 
Finally, U and V in Eq. (JT5J) are the matrices that diagonalize chargino mass matrix, 
which is defined by 

U*M X+ V~ l = diag (m~+, m~+) , (17) 



with 



M 2 \/2mw sin /3 
\f2mw cos /3 /i 



Here M 2 is the weak gaugino mass, fi is the supersymmetric Higgs mixing term, and 
tan/3 = v 2 /vi is the ratio of the vacuum expectation value (VEV) of the two-Higgs 
doublet. From above functions one can see that 

• The second term in Rt L and R g L are enhanced by the large ratio m Xi /mb ~ 30 — 100 
and therefore provide a large chargino contribution to both CV 7 and C 8g . 

• For a large tan/3, say around 30 to 50, the Yukawa coupling Yj, is also large and 
leads to a further enhancement to the LL terms in CV 7 and C$ g . 
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In order to illustrate clearly the impact of the chargino contributions to B — > X s ^f process, 
it is very useful to present the explicit dependence of the Wilson coefficients C 7l (M w ) 
and C$ g (M w ) on the relevant mass insertions. For gaugino mass M 2 = 200 GeV, averaged 
squark mass rh = 500 GeV, fi = 300 GeV and tan {3 = 20, we obtain numerically 



C^(M W ) 
Cj g (M w ) 

C^(M W ) 

CUM W ) 



0.411(^)31 + 1-869(^)32 + 0.002(^)31 + 0.011(5^ 
0.104(^)31 + 0.475(5« L ) 32 + 0.001(^)31 + 0.004(5", 



°RL) 32, 

1.112(^) 31 + 5.062(5^)32 + 0.009(5^) 31 + 0.042(5^) 32 , 
0.507(5^31 + 2.309(5^)32 + 0.006(5^) 3 i + 0.026(5^) 32 . 



(19) 
(20) 

(21) 

(22) 



It is evident that the MIA parameter (5^)3^ (^lr)^) dominates the chargino contribution 
to C 77i 8 9 (CV 7>8g ). However, one should be careful with this contribution since it is also the 
main contribution to the b —>■ 57, and stringent constraints on (5^ L ) 32 are usually obtained, 
specially with large tan/3. Finally, as expected from Eq. ([To]) , only LL contributions to 
C 7l and C$ g have strong dependence on the value of tan j3. With tan (3 = 40, for instance, 
these contributions are enhanced by a factor 2, while the LR part of and C$ g are 
changing from the previous values by less than 2%. 
Now let us turn to the gluino contributions in the b 
at scale nw are given by 



s transition. Its contribution 



to C 7y and C> 



8«,7T 



9V2G F m 2 



A; 



(6 d LL ) 23 M 3 (x § ) + (5 d LR ) 



7Tl~ 

11 if, 



(23) 



V2G F m 2 



-1 



{Sl L )23[^M 3 {x s )+3M A {x- 9t 



rrig I 1 



+(5 d LR h3-T [Mx s )+3M 2 (x~ g ] 



o 77 



r<9 

^8g 



5a, rr 



9V2G F m 2 



a; 



(Kr)2 3 M 3 (x- 9 ) + (5^) 23 -^M!(x,) 
(**n)aa (hUx- g )+3M 4 (x- g ) 



a s iT 



V2G F m 2 



K 1 



(24) 
(25) 



+ (4l) 23 ^ f^M!(x,) + 3M 2 (x,; 
mi) \6 



(26) 



j g ~ "M^j/ih an d the functions Mi(x) can be found in Ref. 



where x & 

Here in order to understand the impact of the gluino contributions to the branching 
ratio of b — > 57, we also present the explicit dependence of the Wilson coefficients CV 7)i 
on the relevant mass insertions. For m = 500Gev, and x = 1, we obtain 



C 9 7j (M w ) 
Cj g {M w ) 
Cl(M w ) 
C! 9 (M W ) 

It is easy to see that the MIA parameter (5^ R ) 2 3 ( (5^ L ) 2 3 
<g (Pig \ 

7 7 ,8g l 7 7 ,8 9 J- 



-0.049(5^)38 - 17.168(^38, 
-0.130(5^) 23 - 64.379(5^) 23 , 
-0.049(5^) 23 - 17.168(5^ L ) 23 , 
-0.130(5^)28 - 64.379(5^ L ) 23 . 



(27) 
(28) 
(29) 
(30) 



dominates the gluino contri- 



bution to H „ {C g 
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After the inclusion of the chargino and gluino contributions, the total Wilson coeffi- 
cients C 7l (M w ) and C 8g (M w ) can be written as 

C 7y(8g) (M w ) = C^% g) (M w ) + C^ {8g) (M w ) + Cl (8g) (M w ), (31) 
C 7 , {89) (M W ) = C^ [8g) (M w ) + Cl m (M w ). (32) 

Since the heavy SUSY particles have been integrated out at the scale fiw — M w , the QCD 
running of the the Wilson coefficients Ci(M w ) down to the lower energy scale fi b = 0{m b ) 
after including the new physics contributions is the same as in the SM. 



III. THE CONSTRAINTS ON THE MASS INSERTION PARAMETERS 

In this section, we will update the constraints on the mass insertion parameters (5^' d )jj 
by considering the data of the branching ratios for the B — > decay and the semilep- 
tonic B — > XJ + l~ decay 

For B — > Xgj decay, the new world average as given by Heavy Flavor Averaging Group 
(HFAG) is 

Br(B -> X sl ) = (3.3918:1?) x 10 ~ 4 - ( 33 ) 

where the error is generally treated as la error. From this result the following bounds (at 
3a level) are obtained 

2.53 x 1(T 4 < Br(B -> X s7 ) < 4.34 x 10" 4 . (34) 

Here the experimental error at 3a level has been added in quadrature with the theoretical 
error in the SM as given in Eq. fjlOj) . From the bounds in Eq. (j34j) . the magnitude-but not 
the sign- of the Wilson coefficient C 7l (m,h) can be strongly constrained. In the SM, the 
allowed ranges for a real C 77 (m fe ) are found to be 

- 0.360 < C 7l (m b ) < -0.248, (35) 
0.454 < C 7j (m b ) < 0.564, . (36) 

But the second range where C 7l {m b ) is positive in sign is strongly disfavored by the 
measured branching ratios of B — > X s l + l~ decays 

By using M 2 = 200 GeV, m = 500 GeV, /i = 300 GeV, tan/3 = 20 and x- g = 1, we find 
the numerical result for the Wilson coefficient C 7l (m b ) after the inclusion of new physics 
contributions from chargino and gluino penguins, 

C 77 (m b ) = -0.3052 - 0.045(5^)23 + 0.284(5^ L ) 31 + 1. 295((^ L )32 

-17.41(5^)23 - 0.002(5« L )3i + 0.008(5^)32, (37) 

where the first term is the C^ M (mh) at the NLO level, other terms represent the new 
physics contributions considered here. Using the same input parameters, the chiral flipped 
part C 7l (m b ) is found to be 

C 7l (m h ) = -0.045(5^)23 + 0.007(5^)31 + 0.030(5^)32 

-17.41(5^)23 + 0.791(5^) 31 + 3.601(5£ H )32. (38) 

One can see from Eqs. lJBTj) and (|HH|) that 
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(i) for the MIA parameters (<^ L )3i,32, ($t L )23, (8rr)31,32 and {5 RR ) 2 3, there is no real 
constraint 2 can be derived from the data of B — > X s 7 and B —>■ X s l + l~ because of 
the tiny coefficients of these terms. We set these parameters zero in the following 
analysis. 

(ii) the MIA parameters {5 d LR ) 2 3 and (i^J 31,32 ((^lr)3i,32 and {8 RL )23 ) dominate 
the new physics contributions to CV 7 (m;,) (CV 7 (m&)). 

When we take the new physics contributions into account, the branching ratio of 
B — > X s 7 decay can be written as 

Br(B -> X s7 ) = Br^(B -> X^)^*^ ■ J^7(z) [\CrM + V(/i b )| 2 



+A{/i b ) + A(/i 6 ) + C 77 (/i 6 ) 



(39) 



where CV 7 (/ift) includes the SM and the SUSY contributions as shown in eq.(p?7|). while 
the last term refers to the new physics contribution coming from the chiral-ffipped oper- 
ators 77 and 8g . For the functions V(/x&), A(fi b ) and A(/i b ) in Eq. (J3TJJ) . the possible 
new physics contributions to these small high order quantities are also small and will be 
neglected in the following numerical calculations. 

Since the gluino and the chargino contributions are given in terms of the MIA parame- 
ters of the up and down squark sectors, they are, in principle, independent and could have 
constructive or destructive interference between themselves or with the SM contribution. 
If we consider all MIA parameters in Eqs. (J37|) and (|38|) simultaneously, no meaningful 
constraint can be obtained indeed. So we consider the cases of one mass insertion and 
two mass insertion only: 

(i) one mass insertion: the SM contribution and one non-zero MIA parameter are 
taken into account; 

(ii) two mass insertion: the SM contribution and two non-zero MIA parameters are 
taken into account; 



A. One mass insertion: gluino contributions 

We now consider the case of only one mass insertion parameter is non-zero at a time. 
The MIA parameters are in general may be complex and can be written as 

(40) 

We firstly consider the case of (6^)23 7^ 0. The MIA parameter (Sf jB ) 2 3 dominate the 
gluino contribution to the Wilson coefficient C7 7 (mb), and therefore strong constraints on 
its magnitude, as listed in Table ITT| can be derived by considering the data of BriB — > 



2 The term "real constraint" means that the constraint of ( 



< 1. 
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TABLE II: Upper bounds on MIA parameters | (5^^)23 1 and | (5^)23 1 from the data of Br(B —> 
Xgj) for m = 500 GeV, Xg = 0.3, 1, 3 and arg(5L R ) 23 = (a), tt/2 (b) and tx (c), respectively. 



X 


0.3 


1 


3 




(a) 0.0024 


(a) 0.0021 


(a)0.0014 


l&)23l 


(b) 0.0099 


(b)0.0083 


(b)0.0054 




(c) 0.0033 


(c) 0.0030 


(c) 0.0019 




0.0102 


0.0107 


0.0103 



X s "f). In numerical calculations, we assume m = 500 GeV, x g = 0.3, 1, 3, and arg(5LR)23 = 
0, tt/2 and tt, respectively. In Table HH we also show the constraints on the magnitude of 
the parameter (6^)23- 

One can see from Table ILT1 that the upper bounds on | (5^)23 1 are sensitive to the phase 
of this mass insertion. But the phase of the parameter (6^)23 is still free parameter. Since 
the parameter (6^)23 appears in CV 7 (m;,) only, which does not mix with CV 7 (m&). 

B. One mass insertion: Chargino contributions 

Now we consider the chargino contribution as the dominant SUSY effect to b — > 57. 
From Eqs. (|3T|) and Eq. ()38|). one can see easily that the LL and LR sector mass insertions 
give the dominant contribution to the Wilson coefficients CV 7 (mh) and CV 7 (mb). 

From the data of Br(B — > X s ^f) as given in Eq. (|34jl and the requirement of a SM-like 
C 77 (m b ) as indicated by the measured branching ratios of B — > X s l + l~ decay, the upper 
bounds on (££2,) 31,32 and (5^)31,32 as the function of the gaugino mass M 2 are given in 
Table. 11111 In the numerical calculation, we use m = 500 GeV, \i = ±300 GeV, and 
tan j3 = 20. For larger values of tan f3, the upper bounds as listed in Table IIHI will be 
scaled by a factor of 20/ tan/3. 

By taking into account the data of BriB — > X s l + l~), the constraints on the MIA 
parameters become more stronger than those obtained by considering the data of BriB — > 
X s 7) only. In Fig. [T] and El we show the the phase dependence of the bounds on |(^l)3iI 
and |(^ L ) 32 | by assuming M 2 = 200 GeV, m = 500 GeV, fi = 300 GeV and tan/9 = 20. 
Here both the shadow regions of A and B are allowed by the data of Br(B — > X B j), but 
the region A is excluded by considering the additional restrictions coming from the data 
of Br(B ^ X s l + l~). 

Furthermore, because of the SU(2) gauge invariance the soft scalar mass Mq is common 
for the up and down sectors. Therefore, one can get the following relations between the 
up and down type mass insertions 

(SLh = [V^km^Vckm]^ ■ (41) 

For the elements ij = 31, 32, we have 

(*2l)31 = (S u LL ) 3l - X(5 U LL )32 + 0{\ 2 ) , (42) 
(^L)aa = (^x)32 + K^llU + 0(X 2 ) . (43) 
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FIG. 1: Plot of 1(^^)31 1 as a function of its phase 9. Here both A and B regions are allowed 
by the data of Br(B — ► X s ^y), but the region A is excluded by the requirement of a SM-like 
C 7j (m b ). 



0. 7 




n/4 n/2 3n/4 

e 



FIG. 2: Plot of 1(^^)32! as a function of its phase 9. Here both A and B regions are allowed 
by the data of Br(B — ► X s j), but the region A is excluded by the requirement of a SM-like 
C 7j (m b ). 
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TABLE III: Upper bounds on |(<5£ L )3i, 32 | and 1 0^)31,32 1 from the data of Br(B -► X s j) and 
Br{B -> X s l + l~) for m = 500 GeV, /i = 300 GeV, tan/3 = 20 and the phase 6 = (a), vr/2 (b) 
and 7T (c), respectively. 





= 300 


H = -300 




200 


400 


600 


200 


400 


600 


1(^)31 1 


(a) 0.185 

(b) 0.618 

(c) 0.127 


(a) 0.267 

(b) 0.894 

(c) 0.185 


(a) 0.365 

(b) 1.230 

(c) 0.255 


(a) 0.135 

(b) 0.532 

(c) 0.192 


(a) 0.191 

(b) 0.765 

(c) 0.282 


(a) 0.264 

(b) 1.05 

(c) 0.382 


l(^L)32| 


(a) 0.042 

(b) 0.138 

(c) 0.028 


(a) 0.057 

(b) 0.195 

(c) 0.042 


(a) 0.082 

(b) 0.270 

(c) 0.056 


(a) 0.028 

(b) 0.116 

(c) 0.042 


(a) 0.042 

(b) 0.168 

(c) 0.061 


(a) 0.058 

(b) 0.230 

(c) 0.084 


1(^)31 1 


0.201 


0.197 


0.195 


0.204 


0.199 


0.197 


1(^)32 1 


0.044 


0.044 


0.043 


0.045 


0.044 


0.043 



TABLE IV: Upper bounds on |(<^ L )i 3 , 2 3| obtained from Table ITTT1 and Eq. (@2J) and (@3J). 





ji = 300 


H = -300 




200 


400 


600 


200 


400 


600 


Ife)l3| 


(a) 0.176 

(b) 0.588 

(c) 0.121 


(a) 0.255 

(b) 0.851 

(c) 0.176 


(a) 0.347 

(b) 1.17 

(c) 0.243 


(a) 0.129 

(b) 0.507 

(c) 0.183 


(a) 0.182 

(b) 0.728 

(c) 0.269 


(a) 0.251 

(b) 0.999 

(c) 0.364 


1(^)23 1 


(a) 0.083 

(b) 0.274 

(c) 0.056 


(a) 0.116 

(b) 0.391 

(c) 0.083 


(a) 0.162 

(b) 0.540 

(c) 0.112 


(a) 0.058 

(b) 0.233 

(c) 0.084 


(a) 0.084 

(b) 0.336 

(c) 0.123 


(a) 0.116 

(b) 0.461 

(c) 0.168 



Consequently, the upper bounds on 1(5^)31,32 1 can be conveyed to a constraint on 
1(5^)31,32! which equals to | (<5^^) 13,23 1 , due to the hermiticity of (M^ll- From the 
bounds as given in Table ITTT1 and the relation in eqs. (jI2J) and (j4*HJ) . the upper bounds 
on 1(5^)13,23! can be obtained, as listed in Table HV] where the terms proportional to A 2 
or higher powers are neglected. These are the strongest constraints one may obtain on 
I (^ll) 13,23 1? an d therefore it should be taken into account in analyzing the "LL" part of 
the gluino contribution to the h — » s and b — *> d transitions. 

C. Two mass insertion 

Finally we consider the scenario in which two mass insertion parameters are non-zero 
at a time. Because of the interference between different terms, the situation become more 
complicated than the cases of one mass insertions. By assuming M 2 = 200 GeV, m = 500 
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0.030 
0.025 
0.020 

m 

"~5 0.015 

to 

0.010 

0.005 

0.000 
0. 



FIG. 3: Contour plot of Br(B — > X s 7) x 10 4 as a function of |(^l)32| an d K^la^I with a same 
zero phase. 

GeV, Xg — 1 and tan/? = 20, we consider only two typical cases: 
Case-A: (S1 L ) 32 and (5£r)23 m C 7l (rrib) are non-zero, i.e. 

CV 7 K) = "0.3052 - 17.41(^) 2 3 + 1.295(5^)32, 

C 7l (m b ) = 0. (44) 

Case-B: (<5^ L )23 and 5^)32 in C 77 (m; ) ) are non-zero, i.e. 
C 7l (m b ) = -0.3052, 

CV 7 (m fe ) = -17.41(^)23 + 3.601(^)32- (45) 

We firstly consider the case A, where the SM and the chargino contributions are taken 
into account. In Fig. |31 we show the contour plot of BriB — > Xg'j) (in unit of 10~ 4 ) as 
a function of | (5^)32 1 and | (5^)23! with a same zero phase. Where the dotted, dashed 
line shows the lower and upper bounds of the branching ratios as given in Eq. (|34p: 
2.53 x 10" 4 < Br{B -> X s ^) < 4.34 x 10" 4 , while the solid line refers to the central 
value of the data. The region between the dotted and dashed lines are allowed by the 
data of both B — > X s 7 and B — > X s l + l~ decay. The Fig. H] shows the same contour plot 
as Fig. (JHJ), but for a non-zero phase of 9 = tt/2. The region between two solid lines is 
allowed by the data of both B — > X s 7 and B — > X s l + l~ decay. 

We now consider the case B, where the SM and the gluino contributions are taken into 
account. In Fig. we show the contour plot of BriB — > X s ^f) (in unit of 10 -4 ) as a 
function of | (5^)32! and |(^ L ) 23 | with a same zero phase. The region between the dotted 
and solid lines are allowed by the data of both B — > X s 7 and B — > X s l + l~ decay. 
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0.15 



0.30 




FIG. 5: Contour plot of Br(B — ► X s ~f) x 10 4 as a function of | (^Eje)32 1 an d K^l^I with the 
same zero phase . 

IV. SUMMARY 

In this paper, by comparing the theoretical predictions with the corresponding mea- 
sured branching ratios of B — > X s ^ and B — > X s l + l~ decays, we make an update for the 
upper bounds on the mass insertions parameters (S^)ij appeared in the MSSM. 

From the well measured B — > X,7 decays, the strong constraints on the MIA param- 
eters can be obtained. From the latest Belle and BaBar measurements of the inclusive 
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B — > X s l + l~ branching ratios, further constraints can be derived. We found that the 
information from the measured branching ratio of B — > X s l + l~ decay can help us to 
improve the upper bounds on the mass insertions parameters {S\ B ) S j i3 m the MSSM. 

We focus on the possible large new physics contributions to the Wilson coefficients 
CV 7 ,8g and Cjy^g coming from chargino and gluino penguin diagrams. Throughout our 
analysis, the SM contributions are always taken into account. From the numerical calcu- 
lations, we found the following points 

1. In the one mass insertion approximation, strong upper bounds on the mass insertion 
parameters (^ L ) 3132 , (^ll) 13 23 ) 

23 are obtained, as 

collected in Tables ITIllTVl 

2. As shown explicitly in Figs. Q and 121 the region A allowed by the data of Br(B — > 
Xs^) is excluded by the requirement of a SM-like CV 7 (m{,) imposed by the data of 
Br[B -> X s l+l~). 

3. Under two mass insertion approximation, strong upper bounds on (^ L ) 23 , and 
{^lr) 23> and (^lr)23 ano - (3rl) 23 can a ^ so De obtained, as shown in Figs. EHS1 
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